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ABSTRACT 



Context. Double-lined, detached eclipsing binaries are our main source for accurate stellar masses and radii. This paper 
is the first in a series with focus on the upper half of the main-sequence band and tests of 1-2 Mq evolutionary models. 
Aims. We aim to determine absolute dimensions and abundances for the detached eclipsing binary V1130Tau, and to 
perform a detailed comparison with results from recent stellar evolutionary models. 

Methods, uvby light curves and uvhyjS standard photometry have been obtained with the Stromgren Automatic 
. Telescope, and high-resolution spectra have been acquired at the FEROS spectrograph; both are ESO, La Silla fa- 

^U,' cilities. We have applied the Wilson-Devinney model for the photometric analysis, spectroscopic elements are based 

on radial velocities measured via broadening functions, and [Fe/H] abundances have been determined from synthetic 
spectra and uvby calibrations. 

Results. V1130Tau is a bright {mv = 6.56), nearby (71 ± 2 pc) detached system with a circular orbit (P = OfSO). The 



o 

^ ' components are deformed with filling factors above 0.9. Their masses and radii have been established to 0.6-0.7%. We 

' derive a [Fe/H] abundance of —0.25 ± 0.10. The measured rotational velocities, 92.4 ± 1.1 (primary) and 104.7 ± 2.7 

(secondary) kms~^, are in fair agreement with synchronization. The larger 1.39 Mq secondary component has evolved 
to the middle of the main-sequence band and is slightly cooler than the 1.31 Mq primary. Yonsai-Yale, BaSTI, and 
^ ' Granada evolutionary models for the observed metal abundance and a 'normal' He content of y = 0.25—0.26, marginally 

QQ ' reproduce the components at ages between 1.8 and 2.1 Gyr. All such models are, however, systematically about 200 

I K hotter than observed and predict ages for the more massive component, which are systematically higher than for 

the less massive component. These trends can not be removed by adjusting the amount of core overshoot or envelope 
• convection level, or by including rotation in the model calculations. They may be due to proximity effects in V1130 Tau, 

• ' but on the other hand, we find excellent agreement for 2.5-2.8 Gyr Granada models with a slightly lower Y of 0.23-0.24. 

I Conclusions. VI 130 Tau is a valuable addition to the very few well-studied 1-2 Mq binaries with component(s) in the 

, upper half of the main-sequence band, or beyond. The stars are not evolved enough to provide new information on the 

dependence of core overshoot on mass (and abundance), but might - together with a larger sample of well-detached 
^— ^ ' systems - be useful for further tuning of the helium enrichment law. Analyses of such systems are in progress. 
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1. Introduction weak. Finally, it is situated at a distance of only 71 pc, 

^ , . 1 n 1 -1 1 1 <• 1 meaning that it belongs to the (small) group of eclipsing 

In^ this paper, we Present the first detailed study of the ^-^^^-^^ ^-^^-^ ^25 pc, discussed by Popper which 

bright my = 6.56) early F-type, double-lmed eclipsing ^ ^ j ^ improving the radiative flux scale. 

binary V1130 Tau. The orbital period is short (P = 0^80 , t ^.u f n • ^ 4. ■ u w • ^ 

, , , , , . . .11 1 , 1 , , r 1 • ; In the loUowmg, we determine absolute dimensions and 

but the system is still detached, and tor several reasons it , , u j i r i, r i ^ j 

. , ,. . •; • 1 1 ;i . abundances, based on analyses 01 new uuDw light curves ana 

IS an interesting case, i'lrst, it is more evolved than most , . , , '. , , TmonT^ j. i 

„ , , 11 i 1- 1 1 n i • i high-resolution spectra, and compare V1130 iau to several 

or the well studied early F-type mam sequence systems; , , w j i o^i i. j. j.u 

^ 11 ii .•' 1 ^ _L 1 ^ modern stellar evolutionary models, ihroughout the paper, 
actually the more massive, larger component has become , ,. j . .i r j t 1 
T 1 , 1 , AT J. -i. • 1. J ^ -u i 1 the component eclipsed at the slightly deeper echpse at 
the slightly cooler one. Mext, it is reported to be metal- , „f.. „ ^, ■ / \ t .t 
° phase 0.0 is referred to as the primary (p), and the other 

Send offprint requests to: J.V. Clausen, as the secondary (s) component. 

e-mail: jvc@nbi.ku.dk 

* Based on observations carried out at the Stromgren 

Automatic Telescope (SAT) and the 1.5m telescope at ESO, La 2. V1130 Tail 

Silla (62.H-0319, 62.L-0284, 63.H-0080, 64.L-0031, 66.D-0178) 

** Table 11 is available in electronic form at the CDS via HD 24133 (CSV 356, H IP 17988) was confirmed to be 

anonymous ftp to cdsarc. u-strasbg.fr (130.79.128.5) or via variable by Olsen (|1983p . Based on the uvby photometry. 



http://cdsweb.u-strasbg.fr/cgi-bin/qcat? J/A-|-A/( Olsen supplied Abt with a list containing about 800 po- 
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Table 1. Photometric data for VllSOTau and the comparison stars. 



Object 


Sp. Type 


Ref. 


V 


(T 


h — y 


a 


mi 


o 


Cl 


(7 


'N{uvby) 




a 


N(/3) 


VllSOTau 


F2 V 


CIO 


6.556 


9 


0.263 


3 


0.140 


7 


0.478 


9 


19 


2.653 


6 


32 






F88 


6.594 


30 


0.272 





0.136 


2 


0.478 


1 


2 


2.652 




1 






083 


6.639 


50 


0.276 


1 


0.124 


4 


0.474 


3 


4 


2.652 


7 


3 


HD23503 


F2/3V=' 


CIO 


8.262 


6 


0.261 


5 


0.174 


8 


0.509 


9 


156 


2.686 


8 


22 




083 


8.251 


1 


0.272 





0.163 


9 


0.510 


1 


2 












094 


8.264 


5 


0.266 


3 


0.167 


4 


0.497 


6 


1 


2.666 


6 


1 


HD24552 


Gl V^' 


CIO 


7.979 


6 


0.386 


5 


0.200 


8 


0.325 


8 


107 


2.596 


8 


16 






083 


7.972 


6 


0.392 


3 


0.200 


1 


0.329 


3 


2 












094 


7.975 


5 


0.387 


3 


0.196 


4 


0.326 


6 


1 


2.575 


6 


1 


HD25059 




CIO 


9.161 


5 


0.391 


4 


0.192 


8 


0.329 


9 


101 


2.591 


9 


19 






093 






0.399 


3 


0.164 


5 


0.359 


6 


1 












094 


9.163 


4 


0.392 


3 


0.185 


5 


0.335 


7 


1 


2.583 


6 


1 



^ Houk & Swift lfT999)) : for VllSOTau, see also spectral types in Sect. [2] 

NOTE 1: References are: CIO = This paper, F88 = Franco (|T988ll . 083 = Olsen (tT983)) . 093 = Olsen (fT993ll . 094 = uvhy. Olsen 
lfT994l) . Olsen (unpublished). 

NOTE 2: For VllSOTau, the uvby information by C09 is the mean value at phases 0.25 and 0.75, and the (3 information is the 
mean value outside eclipses. 

NOTE 3: N is the total number of observations used to form the mean values, and a is the rms error (per observation) in mmag. 



tentially weak-lined A5-G0 stars, and Abt (|1986p subse- 
quently classified HD 24133 as FSVwl (A5 met) and also 
reported it to be a double-lined spectroscopic binary. Gray 
(Pgg)) classified HD 24133 as F5Vm-2, Gray & Garrison 
I1989P confirmed that its metal lines have the strength of 
an A5 star, and Gray et al. (|200ip found it to be a fairly 
rapid rotator, but clearly metal- weak. The eclipsing na- 
ture was discovered by Hipparcos (ESA 119971 orbital pe- 
riod 0'? 7988710), and HD 24133 was subsequently assigned 
the variable name V1130Tau (Kazarovets et al. I1999p . 
Rucinski et al. (|2003p determined a spectroscopic orbit, 
leading to (Mp + Ms)sm^i = 2.41 ± 0.03 Mq, and noted 
that V1130 Tau is one of the shortest period detached early 
F-type systems. Besides a few times of minima, nothing has 
been published on this binary since then. 



3. Photometry 

Below, we present the new photometric material for 
V1130Tau and refer to Clausen et al. (120011 hereafter 
CHOOl) for further details on observation and reduction 
procedures, and determination of times of minima. 

3.1. Light curves for VllSOTau 

The differential uvhy light curves of VI 130 Tau were ob- 
served at the Stromgren Automatic Telescope (SAT) at 
ESO, La Silla and its 6-channel uvby (5 photometer on 
59 nights between October 1997 and November 1998 
(JD2450727-2451120). They contain 583 points per band 
with all phases covered at least twice. The observations 
were done through an 18 arcsec diameter circular di- 
aphragm at airmasses between 1.2 and 1.8. HD 23503, 
HD 24552, and HD 25059 - all within a few degrees from 
V1130Tau on the sky - were used as comparison stars 
and were all found to be constant within a few mmag; 
see Table [T] The light curves were calculated relative to 
HD 23503, but all comparison star observations were used, 
shifting them first to the same light level. The average ac- 



curacy per point is about 4-5 mmag {ybv) and 7 mmag 
(«). 

As seen from Fig. [T] VllSOTau is detached but fairly 
close, with y eclipse depths of about 0.4 mag. Primary 
eclipse is only marginally deeper than secondary, meaning 
that the surface fluxes of the components are nearly iden- 
tical. The light curves (Table 11) will only be available in 
electronic form. 

3.2. Standard photometry for VllSOTau 

Standard uvby (3 indices for V1130Tau (between eclipses) 
and the three comparison stars, observed and derived as 
described by CHOOl, are presented in Table [H The indices 
are based on many observations and their precision is high. 
For comparison, we have included published photometry 
from other sources. In general, the agreement is good, but 
individual differences larger than the quoted errors occur. 

3.3. Times of minima and ephemeris for V1130Tau 

Three times of each of primary and secondary minimum 
have been established from the uvhy light curve observa- 
tions; see Table [21 A list of earlier times of minima was 
kindly provided by Kreiner; see Kreiner et al. (2001) and 
Kreiner (|2004|P ]. Except for two unpublished times based 
on Hipparcos photometry, which showed large deviations, 
they were included in the ephemeris analysis together with 
the recently published time of primary minimum by Brat 
et al. p008|) . 

Assuming a circular orbit, we derive the following linear 
ephemeris from a weighted least squares fit to all accepted 
times of minima: 

Min 1= 2450770.69601 + 0^798868143 x £: . . 

±9 ±38 ^ > 

Separate weighted linear least squares fits to the times of 
primary and secondary minima lead to identical orbital pe- 

^ [http : //www . as . ap . krakow . pi/ ephem) 
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Fig. 1. y light curve and b ~ y and u ~ b colour curves (instrumental system) for VllSOTau. 



Table 2. Times of primary (P) and secondary (S) minima 
of V1130Tau determined from the uvby observations. 



HJD 

- 2 400 000 


rms 


Type 


0-C2 


50742.73603 


0.00040 


P 


0.00041 


50770.69594 


0.00020 


P 


-0.00007 


50778.68480 


0.00020 


P 


0.00011 


50774.73329 


0.00040 


S 


0.00049 


50776.68754 


0.00020 


s 


0.00002 


50780.68191 


0.00020 


s 


0.00005 



" Calculated from the ephemeris given in Eq. [T] 



riods, and we adopt Eq. [l]for the analyses of the uvby light 
curves and radial velocities. 

3.4. Photometric elements 

Since the relative radii of the components of VllSOTau 
are fairly large, 0.25-0.30, we have adopted the Wilson- 
Devinney model (Wilson & Devinney 1971; Wilson 119791 
[T9901ITM31 Van Hamme & Wilson 2003 ) for the hght curve 
analyses. We have used the JKTWD code developed by 



J. SouthworthQ, which is based on the 2003 version of 
the 'Binary Star Observables Program' by Wilson et alH 
This code was recently applied for the light curve analyses 
of DW Car (Southworth & Clausen [20D7|) and V380 Cyg 
(Pavlovski et al. 2009J. 

Mode 2 (detached binaries) was used throughout, and 
the stellar atmosphere approximation functions for the uvby 
bands were adopted (Van Hamme & Wilson [2003( Kurucz 
I1993p . The effective temperature of the primary compo- 
nent was kept at 6500 K; see Sect. |6l A linear limb dark- 
ening law was assumed with coefficients adopted from Van 
Hamme UpM^ . The linear coefficients by Claret CIUUU) are 
about 0.1 larger and lead to a 0?2 lower orbital inclination, 
whereas the radii are practically unchanged. Within errors, 
non-linear limb darkening lead to identical photometric ele- 
ments. Gravity darkening exponents corresponding to con- 
vective atmospheres were applied, and bolometric reflection 
albedo coefficients of 0.5 were chosen, again due to convec- 
tion. The simple reflection mode (MREF — 1) was used; 
we note that the detailed mode gives nearly identical ele- 
ments. The mass ratio between the components was kept 
at the spectroscopic value {Ms/Mp = 1.066 ± 0.004), and 
synchronous rotation was assumed. The light curves were 
analysed independently with at least 10 differential param- 

|http : //www . astro . keele . ac.uk/ ^ jkt7| 
[ftp : / /ftp . astro .ufl . edu/pub/wilson7| 
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Table 3. Photometric solutions for VllSOTau. 





y 


h 


V 


u 


7 (°) 


73 97 


73 84 


73 73 


74 20 


±2 


±4 


±4 


±6 




4.5138 


4.5079 


4.4954 


4.5642 




±39 


±40 


±41 


±64 




4.0727 


4.0753 


4.0775 


4.0805 




±25 


±29 


±30 


±45 


Tp 


0.2956 


0.2961 


0.2973 


0.2911 






u.oooy 


U.oOoD 


U.OOOZ 


k 


1.198 


1.195 


1.189 


1.213 


Tp + Ts 


0.6498 


0.6500 


0.6509 


0.6443 


Up = Us 


0.52 


0.61 


0.67 


0.63 


TcS.3 


6638 


6639 


6643 


6628 


±4 


±3 


±3 


±4 


Ls 1 Lp 


1.444 


1.436 


1.429 


1.436 


a (mmag.) 


4.5 


4.3 


4.3 


6.5 



Note 1: Limb darkening coefficients by van Hamme (1993), 
gravity darkening exponents of 0.33, and bolometric albedo co- 
efficients of 0.5 were adopted, as appropriate for convective en- 
velopes. Tefi,p was assumed to be 6650 K, see Sect. (6] 
Note 2: The errors quoted for the free parameters are the 
formal standard errors determined from the iterative least 
squares solution procedure. 

Table 4. Adopted photometric elements for VllSOTau. 



i 


73? 82 ± 0?20 


Tp 


0.2952 ± 0.0018 


rp{pole) 


0.2862 


rp{point) 


0.3125 


rp{side) 


0.2935 


rp[hack) 


0.3048 




0.3535 ± 0.0021 


Vs (pole) 


0.3369 


Ta (point) 


0.3965 


ra{side) 


0.3506 


rsiback) 


0.3717 



Ls/Lp 1.443 
±23 



1.437 

±23 



1.432 

±23 



1.434 

±23 



Note: The individual luminosity ratios are based on the mean 
stellar and orbital parameters. 



eter corrections, and continuing until they were below 20% 
of the corresponding formal standard errors for all param- 
eters. The stability of the adopted solutions was tested by 
adding 100 more iterations. 

In tables and text, we use the following symbols: i or- 
bital inclination; r relative volume radius; k — rs/vp] U, 
surface potential; u linear limb darkening coefficient; L lu- 
minosity; Teff effective temperature. 

The individual solutions are presented in Table [31 and 
O — C residuals of the h observations from the theoretical 
light curve are shown in Fig. [S] Equally good fits are ob- 
tained in the three other bands, and the rms of the resid- 
uals correspond closely to the observational accuracy. As 
expected, the less massive, slightly hotter primary compo- 
nent is also the smaller one. VllSOTau is detached, but 
the filling factors of the components are above 0.9, and the 
deformation of the secondary component is significant. The 
relative volume radii obtained from the four bands agree 



0,02 
0.01 



Primary eclipse 









0,00 
Phase 



o -0-01 

-0.02 



Secondary eclipse 
— ^ ^ y ^ / ++ ^1 i++> 




-1- ++-('"+ 4_ ++ + + + ^-^ 
+ + + , +^ 





0.50 
Phase 



0-00 



o -0,01 - 
-0.02 - 
0.10 



0.25 
Phase 



o -0.01 
-0.02 



0.65 0.70 



0.75 
Phase 



0.80 0.85 



Fig. 2. (O-C) residuals of the VllSOTau 6-band obser- 
vations from the theoretical light curve computed for the 
photometric elements given in Table [31 



well, except perhaps for the less precise u result for the 
primary. We find no evidence of third light, neither in the 
spectra nor from the light curve solutions, and the small dif- 
ferences in orbital inclinations are probably due to model 
and/or limb darkening effects. The adopted photometric 
elements are listed in Table [H with realistic uncertainties, 
which refiect the formal standard errors and the intera- 
greement of the uvhy results, and also take into account 
the consequences of ±0.1 changes of limb darkening coeffi- 
cients and the uncertainly of the mass ratio. As seen, the 
relative volume radii have been established to about 0.6%. 



4. Spectroscopy 

4.1. Spectroscopic observations 

For radial velocity and abundance determinations, we have 
obtained 18 high-resolution spectra with the FEROS fiber 
echelle spectrograph at the ESO 1.52-m telescope at La 
Silla, Chile (Kaufer et al. [T9991 [2000)) . The spectrograph, 
which resides in a temperature-controlled room, covers 
without interruption the spectral region from the Balmer 
jump to 8700 A, at a constant velocity resolution of 2.7 
kms~^ per pixel (A/AA — 48 000). We refer to Clausen 
et al. (l^niM hereafter CTB08) for details on the reduction 
of the spectra, which were observed between January 1999 
and March 2001; an observing log is given in Table [HI 
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Table 5. Log of the FEROS observations of V1130Tau. 



5 



0.004 r 




-300 -200 -100 100 200 300 

Velocity (km/s) 



HJD-2 400 000^ phase t,;^,^^ S/N^ 



51188.6379 


0.1675 


600 


223 


51207.5446 


0.8344 


600 


171 


51207.5662 


0.8615 


600 


170 


51208.5852 


0.1370 


600 


162 


51209.5918 


0.3971 


600 


200 


51211.5630 


0.8645 


600 


84 


51212.6064 


0.1706 


600 


176 


51385.9274 


0.1288 


600 


142 


51386.9327 


0.3872 


600 


170 


51390.8696 


0.3154 


600 


241 


51391.9171 


0.6266 


600 


225 


51392.8746 


0.8252 


720 


237 


51562.5335 


0.1992 


600 


243 


51562.5962 


0.2777 


600 


281 


51562.6289 


0.3187 


600 


244 


51977.4958 


0.6371 


630 


181 


51978.4897 


0.8812 


600 


179 


51981.4927 


0.6402 


600 


166 



" Refers to mid-exposure 

' Exposure time in seconds 
Signal-to-noise ratio measured around 6070 A. At the shorter 
wavelengths used for the radial velocity measurements it is some- 
what lower 

Table 6. Echelle orders and wavelength ranges used for 
the radial velocity measurements of VllSOTau. 



Order 


Ran 


Se (A) 


Order 


Ran 


je (A) 


55 


4020 


- 4090 


53 


4170 


- 4240 


52 


4220 


- 4330 


50 


4395 


- 4500 


49 


4510 


- 4590 


45 


4905 


- 4975 


44 


4975 


- 5090 


43 


5090 


- 5212 



4.2. Radial velocities 

The radial velocities for VllSOTau were measured from 
eight useful orders (4020 - 5210 A) of the 18 FEROS spec- 
tra, see Table El The selection of this limited number of 
orders was based on initial analyses, which showed that sev- 
eral orders give unreliable results and have to be excluded, 
either because too few lines are available, or because they 
contain defects or are difficult to normalise properly. We 
applied the broadening function (BE) formalism (Rucinski 
fTMOl I2002|. [2004|) . using synthetic templates with no ro- 
tational broadening, calculated for Tcs = 6650 K, log(g) 
= 4.0, and [Ee/H] = —0.25. They were produced with the 
hssynth tool (Bruntt, private communication), which ap- 
plies the SYNTH software (Valen ti fc P iskunov .1996 ) and 
modified ATLAS9 models fHeiter Line information 

was taken from the Vienna Atomic Line Database (VALD; 
Kupka et al. I1999p . BE's were then produced for each of 
the selected orders of each spectrum. 

Due to the deformation of the components (Tabled]) and 
to refiection effects, the lineprofiles and thereby the BE's 
can not a priori be expected to be symmetric. The opti- 
mum method for radial velocity determinations is there- 
fore to fit theoretical BE's which take the proximity effects 
into account; for W UMa systems, see e.g. Rucinski et al 



Fig. 3. Broadening function (thick) obtained for the 5090- 
5212 A region and theoretical fit (thin). The EEROS spec- 
trum was taken at phase 0.834, at HJD=2451207.5446. The 
primary component is to the right. 



P993p and references therein. Such BE's were constructed, 
for both components at each phase and order, from syn- 
thetic rotational broadened line profiles calculated via the 
Wilson-Devinney (WD) model fSect. lH^ . The radial veloc- 
ities were then determined by fitting combined theoretical 
BE's to the observed ones, i.e. by shifting and scaling the 
theoretical BE's for each component until the best fit was 
obtained for the combination. 

Since the observed BE's do not show clear asymmetries, 
we have also applied an approach based on simpler theoret- 
ical BE's, which assume that the stars are spherical, rigid 
rotators (e.g. Kaluzny et al. I2006p . Before the radial veloc- 
ity determination, these BE's were convolved with a Gauss 
profile corresponding to the instrumental resolution. In gen- 
eral, very good fits were obtained from this approach; see 
Fig. El 

In the case of VllSOTau, the two BE approaches re- 
sult in radial velocities for all spectral orders and phases, 
which agree within about 1-2 kms~^, and the differences 
do not correlate with orbital phase and/or velocity sepa- 
ration. In Sect. 14.31 we present orbital solutions from both 
sets of radial velocities. 

As described by e.g. Kaluzny et al. (5006), the pro- 
jected rotational velocities wsini of the components and 
(monochromatic) light /luminosity ratios between them can 
also be obtained from analyses of the simple broadening 
functions mentioned above. We have tested this on syn- 
thetic binary spectra with input rotational velocities of 
90.0 (primary) and 110.0 (secondary) kms~^, correspond- 
ing closely to pseudosynchronous rotation, and a light ratio 
of 1.44, and we find that the method is safe for VllSOTau. 
The rotational velocities determined from the BE analyses 
are within 1 kms~^from the input values, and the light ra- 
tio is reproduced to high precision. Analyses of the observed 
VllSOTau spectra yield mean tisin i velocities of 92. 4± 1.1 
(primary) and 104.7±2.7 kms~^(secondary), and the mean 
light ratio, 1.44 ±0.02, is in perfect agreement with the re- 
sults from the light curve analyses (Table |4]). No significant 
wavelength/order dependencies are seen. 



6 J.V. Clausen et al.: VllSOTau 

Table 7. Spectroscopic orbital solution for VllSOTau. 



Parameter 


WD based BF 


Symmetrical BF 


Mean velocities 
Adopted 


Adjusted quantities: 








Kp (kms"^) 


158.f4±0.34 


158.33 ±0.37 


158.29 ±0.34 


K, (kms-i) 


f 48.45 ± 0.42 


148.49 ± 0.27 


148.52 ±0.33 


7p (kms"^) 


-f0.98±0.27 


-11.12 ±0.30 


-11.05 ±0.28 


7s (kms~^) 


-lf.03±0.34 


-11.41 ±0.22 


-11.22 ±0.27 


Adopted quantities: 








P (days) 


0.798868f43 


0.798868143 


0.798868143 


T (HJD-2 400 000;^ 


50770.6960f 


50770.69601 


50770.69601 


e 


0.00 


0.00 


0.00 


Derived quantities: 








Mp sin^ i (Mq) 


f.f55± 0.007 


1.157 ±0.005 


1.157 ±0.006 


Mssin^i (Mq) 


f .230 ± 0.006 


1.234 ± 0.006 


1.233 ±0.006 


asini (Rq) 


4.839 ± 0.009 


4.843 ± 0.007 


4.842 ± 0.007 


Other quantities pertaining to the fit: 








Nobs 


18 


18 


18 


Time span (days) 


793 


793 


793 


ap (kms~^) 


1.15 


1.26 


1.17 


as (kms~^) 


1.42 


0.93 


1.14 



" Time of central primary eclipse 



4.3. Spectroscopic elements 

Spectroscopic orbits have been derived through analyses of 
the radial velocities obtained from each of the two BF anal- 
yses of the eight selected orders. Since the components of 
VllSOTau are quite close and deformed, the observed light 
center velocities deviate somewhat from the center of mass 
velocities, which are used to determine the Keplerian or- 
bital parameters. Before analysing the velocities, we have 
therefore for each order applied phase dependent correc- 
tions as calculated from the Wilson-Devinney code near 
the corresponding wavelength range; see Sect. 13.41 At the 
observed phases they range between about —1.2 and +1.4 
kms^^ for the primary component and between about —2.2 
and +1.3 kms~^ for the secondary component. Order to or- 
der differences are less than 10% of the corrections; using 
average corrections leads in fact to identical orbital solu- 
tions. 

Next, for each observed phase, mean values of the cor- 
rected radial velocities from the eight selected orders were 
formed, and spectroscopic orbits were then calculated using 
the method of Lehman-Filhes implemented in the SBOIO 
program (Etzel [2"004p . which is a modified and expanded 
version of an earlier code by Wolfe, Horak & Storer (1967). 
A circular orbit was assumed, and the period P and epoch 
T were fixed at the ephemeris values (Eq.[T]). Equal weights 
were assigned to the radial velocities, and the two compo- 
nents were analysed independently (SBl solutions). The el- 
ements are listed in the first two columns of Table [71 and as 
seen, the results from the WD based and the simpler sym- 
metrical BF's agree very well, giving minimum masses accu- 
rate to about 0.6%. For both set of velocities, SB2 analyses 
yield identical semiamplitudes. Within errors, the system 
velocities agree, even without accounting for the small dif- 



* Spectroscopic Binary Orbit Program, 

|http : //mintaka. sdsu. edu/f acuity/ etzel/^ 



ference in gravitational redshift for the components, about 
0.06 kms"^ 

As a further check, we have analysed the eight orders in- 
dependently. The individual semiamplitudes differ slightly 
more than their typical mean errors of 0.6 kms^^, but for 
both BF methods, their mean values agree very well with 
the results presented in Table [71 Finally, applying instead 
mean radial velocities weighted according to the quality of 
the individual order solutions, and/or weighting the mean 
radial velocities according to the S /N ratio of the observed 
spectra (Table [S]), lead to practically identical elements. 
Also, shifting first the velocities from each order by the dif- 
ference between its system velocity and the mean system 
velocity (primary and secondary components treated indi- 
vidually) does not change the elements significantly. 

Based on the results mentioned above, we believe that 
the radial velocity differences from the two BF approaches, 
which for the mean values are within ±1 kms^^, are more 
likely due to imperfections in the observed BF's, affect- 
ing the theoretical fits differently, than to measurable (line) 
asymmetries. Furthermore, the quality of the two datasets 
are comparable with about the same order-to-order spread 
of the velocities. We have therefore taken the pragmatic de- 
cision to base the final orbital elements on the mean values 
of their velocities; see Table [XT] These elements are listed 
in the third column of Table [71 and the corresponding orbits 
are shown in Fig.[H Finally, we note that if the light center 
velocities are applied without corrections, both semiampli- 
tudes become 1.1 kms^^ smaller than listed in Table[71 and 
the derived masses become about 0.03 Mq lower. 

Our results differ slightly from those by Rucinski et al. 
(00051) . Kp = 160.11 ± 0.74 kms-i, = 147.21 ± 0.63 
kms~^, and 7 — — 12.74±0.46 kms^^, and are more accu- 
rate. 
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Fig. 4. Spectroscopic orbital solution for VI 130 Tau (solid 
line: primary; dashed line: secondary) and corrected radial 
velocities (filled circles: primary; open circles: secondary). 
The dotted line (upper panel) represents the mean center- 
of-mass velocity of the system. Phase 0.0 corresponds to 
central primary eclipse. 
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Fig. 5. Part of FEROS spectrum taken at phase 0.168 
(gray) and synthetic binary spectra (thin) calculated for 
[Fe/H] = -0.25 (upper) and [Fe/H] = 0.00 (lower). 

5. Chemical abundances 

Due to the high rotational velocities of the components 
(Sect. [4?2| . a detailed chemical analysis of VllSOTau based 
on the FEROS spectra is difhcult, see Fig. El First, the lines 
with intrinsic equivalent widths below about 100 mA, which 
should preferably be used, are shallow and broad and there- 
fore impossible to measure accurately. Next, line blending 
becomes a serious issue, and finally proper normalization of 
the spectra is difficult, especially in the blue spectral region. 



Primary Secondary 



Absolute dimensions: 






M/Mq 


1.306 ± 0.008 


1.392 ± 0.008 


R/Rq 


1.489 ±0.010 


1.782 ±0.011 


logs (cgs) 


4.208 ± 0.006 


4.080 ± 0.006 


vs'mi^ (kms~^) 


92.4 ± 1.1 


104.7 ±2.7 


Vsync^ (kmS~^) 


90.6 ± 0.6 


108.5 ± 0.7 


Photometric data: 






V 


7.526 ± 0.014 


7.128 ±0.011 


ib-y) 


0.260 ± 0.004 


0.265 ±0.004 


mi 


0.141 ±0.008 


0.140 ±0.008 


Cl 


0.481 ±0.010 


0.476 ± 0.010 




0.000 ± 0.008 


TcS 


6650 ± 70 


6625 ± 70 




3.27 ±0.05 


2.89 ±0.05 


log L/Lq 


0.59 ±0.02 


0.74 ±0.02 


BC 


0.02 


0.01 


Mv 


3.25 ±0.05 


2.88 ±0.05 


Vo - Mv 


4.28 ±0.06 


4.25 ±0.06 


Distance (pc) 


71.6 ±2.1 


70.8 ±2.1 


Abundance: 






[Fe/H] 


-0.25 


±0.10 



Observed rotational velocity 
^ Projected equatorial velocity for synchronous rotation 

Note: Bolometric corrections (BC) by Flower (1996 ) have been 
assumed, together with T^jjq = 5780 K, BCq = —0.08, and 
M,,tQ = 4.74. 

Line by line analyses of either the observed spectra or the 
reconstructed component spectra calculated from disentan- 
gled spectra have therefore not been attempted. We refer 
to CTB08 and Clausen et al. ((M?^ for details on line by 
line analyses of binaries. We have instead established upper 
and lower limits for the metal abundance of VllSOTau by 
comparing the observed spectra and synthetic binary spec- 
tra calculated for a range of scaled solar compositions. The 
synthetic spectra were produced as described in Sect. 14.21 
The overall result, based on inspection of several spectra 
and orders, is that synthetic spectra for metal abundances 
between -0.S5 and -0.15 dex fit the observed spectra equally 
well, whereas e.g. the lines/lineblends for solar abundance 
spectra, as illustrated in Fig. O are clearly too strong. 

In addition, abundances have been derived from vari- 
ous uvby calibrations and the indices listed in Tables [T] and 
E The Holmberg et al. (i2007|) calibration gives [Fe/H] = 
—0.25 ± 0.12 for both components, whereas the 'blue' cali- 
bration by Nordstrom et al. (|2004p gives [Fe/H] = -0.S4 ± 
0.14. For comparison, the older calibration by Edvardsson 
et al. (fTOMI) gives [Fe/H] = -0.27 ± 0.11. 

In conclusion, we confirm that VllSOTau is (slightly) 
metal-weak, see Sect. [2l and adopt [Fe/H] = -0.25 ± 0.10. 

6. Absolute dimensions 

Absolute dimensions for the components of VllSOTau are 
calculated from the elements given in Tables d] and [71 As 
seen in Table [H both masses and (volume) radii have been 
established to an accuracy of 0.6-0.7%. 
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Individual standard uvby indices are included in 
Table [51 as calculated from the combined indices of 
VllSOTau outside eclipses (Table[T]) and the luminosity ra- 
tios (Table SI . According to the calibration by Olsen (|1988p 
and the combined uvbyj3 indices at phase 0.25, there is no 
significant interstellar reddening. 

The adopted effective temperatures (6650 K, 6625 K) 
were calculated from the calibration by Holmberg et al. 
(|2007|) . assuming [Fe/H] = -0.25 (Sect. [5]). The uncertain- 
ties include those of the uvby indices, E{b — y), [Fe/H] 
and the calibration itself. Identical temperatures are ob- 
tained from the calibration by Ramirez & Melendez (2005P, 
whereas that by Alonso et al. (1996) leads to 100 K lower 
values. 2MASS photometry at phase 0.79, where V = 6.555, 
and the V — Kg calibration by Masana et al. p006p gives 
an average temperature of 6600 K. 

The measured rotational velocities (wsin?) are close to 
the projected synchronous velocities. We note that for an 
orbital inclination of 'only' 73?8, the true equatorial veloc- 
ities are about 4% higher. The turbulent dissipation and 
radiative damping formalism of Zahn (19Z7j 1989) predicts 
synchronization times scales of 8.7 x 10^ yr (primary) and 
3.1 X 10^ yrO.l Gyr (secondary), and a time scale for circu- 
larization of 1.2 x 10^ yr. 

The distance to VllSOTau was calculated from the 
'classical' relation (see e.g. CTB08), adopting the solar val- 
ues and bolometric corrections given in Table [S] and ac- 
counting for all error sources. Other BC scales (e.g. Code 
et al. Bessell et al. [TOMl Girardi et al. I^gfe|) give 

nearly identical results. As seen, the distances obtained for 
the two components agree well. The mean distance, 71.2 
pc, which has been established to about 3%, is close to the 
result from the new Hipparcos reduction by van Leeuwen 
(|2007p . 69.8 ± 2.3 pc, but is marginally larger than the orig- 
inal Hipparcos result 65.2 ± 3.3 pc (ESA [T997|) . Finally, we 
note that V1130Tau belongs to the group of eclipsing bina- 
ries within 125 pc, discussed by Popper (|1998p . which could 
be useful for improving the radiative flux scale. 



7. Discussion 

Below, we first compare the absolute dimensions obtained 
for V1130Tau with properties of recent theoretical stellar 
evolutionary models, and we then discuss V1130Tau to- 
gether with the few other similar well-studied eclipsing bi- 
naries available. 

7.1. Comparison with stellar models 

Figs.iniEl andOillustrate the results from comparisons with 
the Yonsei-Yale (F^) evolutionary tracks and isochrones by 
Demarque et al. (20040 They include core overshoot where 
Aos = ^ov/Hp depends on mass and also takes into ac- 
count the composition dependence of M^°™ 0. The mixing 
length parameter in convective envelopes is calibrated using 
the Sun, and is held fixed at I /Hp = 1.7432. The enrich- 
ment law Y — 0.23 -f- 2Z is adopted, together with the solar 
mixture by Grevesse et al. (|1996p . leading to [X,Y,Z)q = 



http : / /www . astro . yale . edu/demarque/yystar . html 

® Defined as "the mass above which stars continue to have 
a substantial convective core even after the end of the pre-MS 
phase." 




7000 6500 6000 

Te,, [K] 

Fig. 6. VllSOTau compared to models calculated for 
[Fe/H] = —0.25. Tracks for the component masses (full 
drawn, thick) and isochrones for 0.5-3.0 Gyr (dashed, step 
0.5 Gyr) are shown. The uncertainty in the location of 
the tracks coming from the mass errors are indicated (full 
drawn, thin). To illustrate the effect of the abundance un- 
certainty, tracks (dotted, blue) and the 2.2 Gyr isochrone 
(dash-dot, blue) for [Fe/H] = —0.15 arc included. 




1.28 1.30 1.32 1.34 1.36 1.38 1.40 



Fig. 7. VllSOTau compared to Y^ models calculated for 
[Fe/H]= -0.25. Isochrones (full drawn) for 1.5, 2.0, 2.13, 
and 2.5 Gyr are shown. To illustrate the effect of the abun- 
dance uncertainty, 2.0 Gyr isochrones for [Fe/H] = —0.15 
(dotted) and [Fe/H] = —0.35 (dashed) are included. 



(0.71564,0.26624,0.01812). A brief description of other as- 
pects of their up-to-date input physics in given by GTB08. 
Only models for [a/Fe]= 0.0 have been included in the 
figures. We have used the abundance, mass, and age inter- 
polation routines provided by the Y^ group. 

As seen from Fig. |6l models for the observed masses 
and abundance, [Fe/H] = —0.25, equivalent to {X ,Y ,Z)q 
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Fig. 8. VllSOTau compared to models calculated for 
[Fe/H] = -0.25. Isochrones (full drawn) for 1.5, 2.0, and 
2.5 Gyr are shown. To illustrate the effect of the abundance 
uncertainty, 2.0 Gyr isochrones for [Fe/H] = —0.15 (dotted) 
and [Fe/H] = —0.35 (dashed) are included. 




1.28 1.30 1.32 1.34 1.36 1.38 1.40 
M/Mq 



Fig. 9. VllSOTau compared to the models listed in 
Table [S] Isochrones for the average ages inferred from 
masses and radii are shown. Y^: thick full drawn (black). 
BaSTI overshoot: thin full drawn (blue). BaSTI standard: 
dotted (blue). Victoria- Regina: dashed (red). 



= (0.7385,0.2510,0.0105), are about 200 K hotter than 
observed. The uncertainty of [Fe/H] is ±0.10 dex, and 
tracks for [Fe/H] = -0.15, equivalent to {X,Y,Z) = 
(0.7310,0.2560,0.0130), fit the components at an age of 
about 2.2 Gyr. This can also be reached for [Fe/H] = —0.25, 
if a slight hypothetical a-element enrichment of [a/Fe] — 
0.15 is introduced. The more massive secondary component 
has evolved to the middle of the main sequence band. 

From a binary perspective, the most fundamental com- 
parison is that based on the scale-independent masses 



Table 9. Models information and average ages inferred 
from masses and radii; see Fig. [51 



Grid 


[Fc/H] 


Y 


Z 


Age (Gyr) 


Yonsai-Yale (F^) 


-0.25 


0.2510 


0.0105 


2.13 


Victoria- Regina 


-0.29 


0.2574 


0.0100 


1.98 


BaSTI (overshoot) 


-0.25 


0.2590 


0.0100 


1.84 


BaSTI (standard) 


-0.25 


0.2590 


0.0100 


1.80 



and radii, as shown in Fig. [T) The [Fe/H] = —0.25 model 
isochrone for 2.13 Gyr marginally fits both components, but 
within the abundance uncertainty, the general trend is that 
the Y^ isochrones predict a higher age for the secondary 
component than for the primary. Although less evident, this 
is also seen in the mass-luminosity diagram (Fig. [5]). 

In Fig. [S] we have included mass-radius comparisons 
with the Victoria- Regina (VRSS grid; VandenBerg et al., 
[5006j!|and BaSTI (Pietrinferni et al., 2004j!| models, which 
differ from Y^, e.g. with respect to input physics. He enrich- 
ment law, and core overshoot treatment. We refer to CTB08 
for a brief description. Basic parameters for the models, all 
with solar scaled abundances, are given in Table [H Like 
the Y^ models, both the standard and overshoot BaSTI 
models marginally fit both components, but at a lower age. 
However, the Victoria- Regina models do not fit VllSOTau 
well. To us, this is surprising, because these models are care- 
fully calibrated by cluster and binary observations. Models 
with [a/Fe] = 0.3 (VR2A grid) can reproduce VllSOTau at 
an age of about 2.15 Gyr, but only for [Fe/H] around -0.40 
dex. 

Thus, except for the Victoria-Regina models, all the 
models with solar scaled abundances we have tested are 
marginally able to reproduce VllSOTau, but we see two 
general trends: First, models for the observed [Fe/H] are 
about 200 K too hot. Second, they systematically predict 
higher ages for the more massive secondary component than 
for the primary. In order to look in more detail into this, 
we have calculated dedicated models for the component 
masses with various parameters tuned. For all models, we 
have adopted Z = 0.010, which is equivalent to the ob- 
served [Fe/H]. We have apphed the Granada code by Claret 
(|2004p . which assumes an enrichment law of y = 0.2A+2.0Z 
together with the solar mixture by Grevesse & Sauval 
(ITM5D . leading to {X,Y,Z)q = (0.704,0.279,0.017). The 
envelope mixing length parameter needed to reproduce the 
Sun is l/Hp = 1.68. The amount of core overshooting is 
given, in units of the pressure scale height, by aov 

Table [TU] lists the models we have investigated. As seen 
in Fig. [TOl the overshoot models (1) with Y calculated 
from the adopted enrichment law are too hot, as seen for 
the other grids. The same is true for models without over- 
shoot (2, not shown). Since the components of VllSOTau 
are rotating quite fast (Table [5]), we have calculated mod- 
els which include rotation as described by Claret (199S|). 
Angular velocities for the models were tuned to reproduce 
the observed equatorial rotational velocities of the compo- 
nents. As expected, such models (3) are cooler than similar 
ones without rotation (1), but the effect is small compared 
to the about 200 K discrepancy. Next, the components of 



http : //wwwl . cadc-ccda.hia-iha.nrc-cnrc .gc . ca/cvo, 
community /VictoriaReginaModels/ 
^ [http : //www . te . astro . it/BASTI/index .php| 
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Table 10. Information on the Claret models and ages inferred from radii; see Figs. [TUlandITT] 



Model/ 


Y 


Z 




O^ov 


Rotation 


Age (Gyr) 


Age (Gyr) 


Linestyle 












Primary 


Secondary 


1 thin, blue 


0.260 


0.010 


1.68 


0.20 


NO 


2.08 ±0.07 


2.24 ±0.03 


2 


0.260 


0.010 


1.68 


0.00 


NO 


2.01 ±0.07 


2.14 ±0.03 


3 dotted, blue 


0.260 


0.010 


1.68 


0.20 


YES 


2.01 ±0.07 


2.19 ±0.03 


4 dashed, blue 


0.260 


0.010 


1.50 


0.20 


NO 


1.98 ±0.07 


2.20 ± 0.03 


5 thick, black 


0.240 


0.010 


1.68 


0.20 


NO 


2.51 ±0.07 


2.60 ± 0.04 


6 


0.230 


0.010 


1.68 


0.20 


NO 


2.75 ±0.07 


2.81 ± 0.04 




Fig. 10. VllSOTau compared to Claret models for the 
observed masses and [Fe/H] abundance. See Table [TUl for 
details and linestyles/colours. 

VllSOTau have (thin) outer convection zones, and we have 
therefore investigated the effect of modifying the envelope 
mixing length parameter. 2D radiation hydrodynamic cal- 
culations by Ludwig et al. (|19991 see also Clausen et al. 
I2009P predict parameters, which are about 0.2 lower than 
for the Sun, and we have therefore adopted 1.50. The mod- 
els (4) become cooler, but again the effect is too small. 
Finally, we have calculated models with a He abundance 
slightly lower than Y = 0.26, as given by the enrichment 
law for Z = 0.01. Tracks for Y = 0.24 (5) and 0.2S (6, 
not shown) actually fit VllSOTau well. If we now turn to 
the ages, as determined from the radii. Fig. [TT] shows that 
these models with lower Y also predict practically identical 
ages for the components. In fact, this also holds if lower aov 
values are adopted; models without overshoot place the pri- 
mary component just at the end of the core hydrogen burn- 
ing phase. On the other hand, all the Granada models for 
Y — 0.26 predict higher ages for the secondary component 
than for the primary, as seen for the other grids. 

Before finishing these comparisons and drawing any def- 
inite conclusions about the need to adjust basic physical 
or chemical ingredients of the models, it is worth remem- 
bering that besides being fast rotating, the components of 
VllSOTau are influenced by their mutual gravitational and 
radiative interactions. They cause not only additional defor- 
mation, but also expansion and some heating, and these ef- 
fects are probably somewhat different for the two stars. We 



Fig. 11. VllSOTau compared to Claret models for the 
observed masses and [Fe/H] abundance. See Table [TUl for 
details. The curves illustrate model radii as function of age 
for the components (upper: secondary; lower: primary). The 
horizontal full drawn lines mark the observed radii of the 
components with errors (dotted lines). The vertical lines 
mark the ages predicted for the secondary component. 

will not elaborate further on the possible implications for 
the model comparisons until additional, similar, but more 
detached binaries have been studied. 

7.2. Comparison with otiier binaries 

Binaries like VllSOTau with component(s) that have 
evolved to the upper half of the main sequence band, or 
beyond, may give important information on core overshoot. 
Already 20 years ago, such systems were found to provide 
strong evidence for convective core overshoot in intermedi- 
ate mass (1.5-2.5 Mq) stars (Andersen et al.[T990). From 
a sample of 2-12 Mq systems, Ribas et al. ()2000|) found 
a significantly increasing of the amount of overshoot with 
increasing stellar mass, whereas Claret (|2007p found that 
it is less pronounced and more uncertain. 

From the onset of core convection up to about 1.5 Mq 
there are, however, only a few relevant, well-studied bi- 
naries: (excluding active systems and systems with nearly 
identical components): AlPhe, BKPeg, BWAqr, and GX 
Gem. Andersen et al. (|1988p found that models without 
core overshoot were able to reproduce AlPhe (1.24 + 1.20 
Mq, components above the main sequence) remarkably 
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well for a normal helium abundanc^EI, whereas Clausen 
P99ip found that models including moderate overshoot 
gave better fits for especially the primary components of 
the slightly more massive systems BWAqr (1.49 + 1.39 
Mq) and BKPeg (1.43 + 1.28 A/q). The latter is consis- 
tent with a lower limit of Uov of about 0.18 for GXGem 
(1.49 + 1.47 Mq), as established by Lacy et al. (2008). 

We had hoped and expected, that VllSOTau could fill 
the mass gap between these systems, but as mentioned in 
Sect. 17. ll this is not the case - Claret models with aov from 

to at least 0.2 can reproduce it perfectly well for Y = 
0.23-0.24. In contradiction to this, Tomasella et al. (2008a', 
I2008b[) report determination of aov from the much younger 
systems V505 Per and V570 Per. 

It is, however, still important to try to calibrate core 
overshoot better from its onset to say 2 Mq. For the 
Victoria-Regina model grids, VandenBerg et al. (2006) 
adopt, from binary and cluster information, a mass and 
abundance dependent amount, setting in around 1.1 Mq 
and gradually increasing up to about 1.7 Mq. Demarque et 
al. ()2004) apply a different ramping algorithm for the 
isochrones, as do Pietrinferni et al. (|2004p for the BaSTI 
calculations. These recipes, and others, need further em- 
pirical tests, and we plan to address that issue in forthcom- 
ing re-analyses of BW Aqr and BK Peg, which will include 
abundance determinations, as well as through new com- 
plete analyses of ALLeo, HD76196, and possibly also the 
NGC752 member DS And. 

Another important aspect is the He abundance and the 
helium-to-metal enrichment ratio, and, through extrapola- 
tion, the primordial He/H abundance ratio. As discussed 
in Sect. 17. li VllSOTau points towards a lower He abun- 
dance and/or enrichment ratio than the four different Y, Z 
prescriptions adopted by the model grids studied. We refer 
to Casagrande et al. (200J7j) for at recent determination of 
AY/AZ based on K dwarfs (2.1±0.9), to Blaser (|2W61) for a 
HII based study (1.41 ±0.62) with references to a variety of 
methods and results, and to Ribas et al. (|2000p and Claret 
& Willems (2002) for determinations based on samples of 
eclipsing binaries (2.2 ±0.8 and 1.9 ±0.6, respectively). We 
believe binaries can give an even better constraint, provided 
detailed heavy element abundance determinations become 
available for a sufficiently large sampled- Such investiga- 
tions are in progress for several systems, and we will return 
to this matter in forthcoming papers. 

Here, we close the issue with a brief historical remark: 
The use of binaries to determine the hydrogen content of 
stars was pioneered by Eddington (|1932p and Stromgren 
(|19S21 11933P , and a few years later Stromgren (1938) also 
used binaries in his classical discussion of the helium con- 
tent of the interior of the stars. Later, binary based helium- 
hydrogen abundance ratio determinations (for Population 

1 stars) were published by Stromgren p967p and Popper et 

al. mm . 

8. Summary and conclusions 

From state-of-the-art observations and analyses, precise 
(0.6-0.7%) absolute dimensions have been established for 
the nearby, early F-type, double-lined, detached eclipsing 

^ see also Torres et al. 1^751 

see Torres et al. (12009 1) Table 3 for the limited material avail- 
able today 



binary V1130Tau. From synthetic spectra and uvby cal- 
ibrations, a metal abundance of [Fe/H] = —0.25 ± 0.10 
has been derived. The 1.39 Mq secondary component has 
evolved to the middle of the main-sequence band and is 
shghtly cooler than the 1.31 Mq primary. The P = (fiSO 
period orbit is circular and the observed rotational veloci- 
ties of the components, 92.4±1.1 (primary) and 104. 7± 2.7 
(secondary) kms~^, correspond closely to synchronization. 

Yonsai-Yale, BaSTI, and Granada evolutionary models 
for the observed metal abundance and a 'normal' He con- 
tent of y = 0.25 — 0.26, as established from the adopted 
helium enrichment laws, marginally reproduce the compo- 
nents at ages between 1.8 and 2.1 Gyr. All such models are, 
however, systematically about 200 K hotter than observed, 
and predict ages for the more massive component, which 
are systematically higher than for the less massive com- 
ponent. The latter is even more pronounced for Victoria- 
Regina models. The two trends can not be removed by ad- 
justing the amount of core overshoot or envelope convection 
level, or by including rotation in the model calculations. 
They may be due to proximity effects in V1130Tau, but 
on the other hand, we find excellent agreement for 2.5-2.8 
Gyr Granada models with a slightly lower Y of 0.23-0.24. 

We had expected that V1130Tau is sufficiently evolved 
to provide new information on the level of core overshoot 
in the 1.1-1.7 Mq interval, where it is believed to ramp 
up, but this is not the case. V1130Tau can be reproduced 
by models calculated for aov from 0.0 to at least 0.2. The 
preference for a helium content of 0.23-0.24 is interesting, 
but more well-detached systems with measured metal abun- 
dances are needed for any firm conclusions on the implica- 
tions for example for the helium enrichment law. We will re- 
turn to these issues in forthcoming papers on other systems 
recently observed within the Copenhagen binary project. 
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Table A.l. Mean radial velocities for VllSOTau. 
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Mean values of the measured radial velocities (kms ^) 
^ Mean correstions. Approximate center of mass velocities are obtained by subtracting these corrections (kms^^) from the 
measured velocities. For the orbital solutions, individual corrections were used for each order. 
O — C residuals (kms~^) from the adopted spectroscopic orbits (Table[7|. 



